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ABSTRACT 

We use redshift observations of two deep 1.4 GHz fields to probe the evolution of the bright end of the 
radio galaxy luminosity function to z — 1.5. We show that the number of galaxies with radio power that 
would correspond locally to an ultraluminous infrared galaxy (ULIG) evolves rapidly over this redshift 
range. The optical spectra and X-ray luminosities are used to refine this result by separating the sources 
with clear active galactic nucleus (AGN) signatures from those that may be dominated by star formation. 
Both populations show extremely rapid evolution over this redshift range. We find that the number of 
sources with ULIG radio power and no clear AGN signatures evolves as (1 + z)^. 

Subject headings: cosmology: observations — galaxies: distances and redshifts — galaxies: active — 
X-rays: galaxies — galaxies: formation — galaxies: evolution 
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1. INTRODUCTION 

Star-forming ultraluminous infrared galaxies, or ULIGs 
{LpiR > 10^^ Lq from Sanders & Mirabel 1996; here we 
mean galaxies without clear signatures of active galactic 
nucleus [AGN] activity in their spectra, even if there may 
be some AGN contribution to their luminosities) , represent 
the most luminous tip of the star-forming galaxy popu- 
lation. Locally, ULIGs are very scarce (Kim & Sanders 
1998), but if, as is widely believed, submillimeter galax- 
ies are the high redshift analogs, then by z = 1 they are 
the dominant contributors to the star formation history 
(e.g., Lilly et al. 1999; Barger, Cowie, & Richards 2000; 
Gispert, Lagache, & Puget 2000; Chapman et al. 2003). 
It is extremely hard to trace the evolution of the ULIG 
population at low redshifts [z < 1.5) because of the diffi- 
culties of mapping large areas with current submillimeter 
detectors and of identifying the optical counterparts to 
low spatial resolution submillimeter or far-infrared (FIR) 
observations. 

Ultradeep decimetric radio surveys can play a major, 
complementary role. The well-known tight correlation 
between FIR luminosity and radio power in local star- 
forming galaxies and radio-quiet AGNs (Condon 1992) 
means that we can identify ULIGs out to substantial red- 
shifts based on their radio power. For a 1.4 GHz sensitivity 
limit of 40 /xJy, an Lpm = 10^^ Lq source will be detected 
to z ~ 1.5. Thus, sources with radio power that would cor- 
respond locally to a ULIG may be seen directly to large 
redshifts. Since the ultradeep 1.4 GHz samples cover large 
fields, do not suffer from extinction, and provide subarcsec- 
ond positional accuracy, it is easy to develop large samples 
with highly complete optical counterpart identifications. 



By contrast, optically selected samples may omit dusty 
sources at these redshifts, while X-ray samples will select 
AGN dominated sources. 

Current decimetric surveys do not probe the normal 
star- forming galaxy populations at z = 1. To the deep- 
est 1.4 GHz sensitivity limits of ^ 20 /iJy (5(t), a source 
with radio power corresponding to Lpm = 10^° Lq would 
only be seen to z ~ 0.3. Thus, attempts to provide maps 
of the entire star formation history over the z = — 1 
range based on decimetric observations (e.g.. Cram 1998; 
Mobasher et al. 1999; Haarsma et al. 2000, hereafter HOO) 
have to rely on model luminosity functions to extrapo- 
late the measured data and obtain star formation rates. 
In these descriptions, much of the star formation lies in 
sources that arc not directly measured by the observations 
(e.g., HOO). 

We use ultradeep 1.4 GHz observations of the Hub- 
ble Deep Field-North (HDF-N; Richards 2000) and SSA13 
(Fomalont et al. 2004) fields to analyze the evolution of the 
radio luminosity function (LF) over the range z = — 1.5. 
Our sample of 346 radio sources represents an increase of 
nearly a factor of five over HOO and enables us to provide 
a much more accurate determination of the z = 1 radio 
LF. We use optical spectral classifications and X-ray char- 
acteristics to separate AGN dominated sources from those 
that may be star formation dominated and to estimate the 
number density of ULIGs. An analysis of the submillime- 
ter and far infrared properties of the sample will appear in 
Wang, Barger, and Cowie 2004 in preparation. We assume 
= i, f^A = |, and i/o = 65 km s^^ Mpc^^. 



2. SAMPLE 



^ Based in part on data obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the the California 
Institute of Technology, the University of California, and NASA and was made possible by the generous financial support of the W. M. Keck 
Foundation. 

^Based in part on data obtained at the VLA, which is a facility of the National Radio Astronomy Observatory (NRAO). The NRAO is a 
facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. 

^Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822. 
^Department of Astronomy, University of Wisconsin-Madison, 475 North Charter Street, Madison, WI 53706. 
^Department of Physics and Astronomy, University of Hawaii, 2505 Correa Road, Honolulu, HI 96822. 
^National Radio Astronomy Observatory, Charlottesville, VA 22903. 
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Richards (2000) presented a catalog of 1.4 GHz sources 
detected in a Very Large Array (VLA) map centered on 
the HDF-N that covers a 40' diameter region with an 
effective resolution of 1.8". The 5a completeness limit 
for compact soiirces in the central region of the map is 
40 iiJy. Fomalont et al. (2004) presented a similar catalog 
within a 17' radius of the 1.4 GHz VLA map {5a com- 
pleteness of 25 fiJy) of the SSA13 field (R.A.= 13'*23'", 
decl.= 42°38"). In both cases, the completeness limits 
for more extended sources arc several times higher. The 
absolute radio positions arc known to 0.1" — 0.2" rms. 
Optical imaging data for the two fields can be found in 
Capak et al. (2003) and Barger et al. (2001), respec- 
tively. Alexander et al. (2003) presented 2 Ms Chandra 
Deep Field-North (CDF-N) X-ray images that cover about 
460 arcmin^ around the HDF-N. Near the aim point, the 
data reach limiting fluxes of « 1.5 x 10"^'^ (0.5 - 2 keV) 
and « 10"^*^ ergs cm~2 s~^ (2-8 keV). Mushotzky et al. 
(2000) and Barger et al. (2001) presented 100 ks Chandra 
exposures of the SSA13 field with limiting fluxes 20 times 
higher than the CDF-N. 

For the HDF-N region, we use a 9.3' radius area 
around the approximate Chandra image center R.A.= 
12'* 36" 51.20" and decl.= 62°13' 43.0" (J2000.0). The 
9.3' radius is the largest circle that can be drawn within 
the Chandra area. 195 of the 372 radio sources lie within 
this circle. For the SSA13 field, we use a rectangular area 
of 175 arcmin^ that contains 221 sources. We have so 
far spectroscopically observed 154 of these sources, cho- 
sen without regard to their radio or optical properties; 
thus, the effective area is 120 arcmin^. Together these two 
restricted samples constitute our radio sample with 346 
sources. All lie well within the half power radius of the 
radio observations, so the two fields should have relatively 
uniform flux selections. 
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Fig. 1. — Redshift versus R magnitude for our radio sample, 
excluding the two spectroscopically identified stars {solid squares — 
spectroscopic redshifts; open triangles — photometric redshifts; large 
solid squares — broad- line AGNs) . Dashed horizontal line shows the 
redshift {z = 1.05) for which a source with the median radio flux of 
80 /iJy would have a radio power that would correspond locally to 
a ULIG. Solid histogram shows the spectroscopic completeness: at 
R < 23, 96% (152 of the 162 sources) are identified; at it = 23.5-24, 
just over 50%; at fainter magnitudes, only a handful with strong 
emission lines. Dashed histogram shows the completeness of the 
combined photometric and spectroscopic redshifts in the HDF-N: 
92% of the il < 25 sources have either a spectroscopic or a photo- 
metric redshift. 



Wc observed our radio sample with the Deep Extra- 
galactic Imaging Multi-Object Spectrograph (DEIAIOS; 
Faber et al. 2003) on the Keck II 10 m telescope and the 
Low-Resolution Imaging Spectrograph (LRIS; Oke et al. 
1995) on Keck I. We spectroscopically identified 192 of 
the 346 sources (55%), two of which are stars. Figure 1 
shows spectroscopic redshift versus R magnitude {solid 
symbols) for our radio sample, as well as the complete- 
ness of our spectroscopic identifications versus R magni- 
tude {solid histogram). A source with the median radio 
flux of 80 /iJy has a radio power corresponding to a lo- 
cal ULIG at z = 1.05 {dashed horizontal line); this cor- 
responds to R magnitudes in the range 22.5 to 25 (see 
Fig. 1). Spectroscopic identifications can be obtained for 
the majority of sources to i? = 24 and for some in the 
i? = 24 to i? = 25 range (see solid histogram) . 

We can also estimate redshifts from photometric colors. 
For the HDF-N, we use the photometric redshifts of Ca- 
pak et al. (2004), which are based on ultradeep imaging in 
seven colors from U to HK' (Capak et al. 2003) and were 
computed using the Bayesian code of Benitez (2000) . This 
code assigns odds that the redshift lies close to the esti- 
mated value; we only use cases with odds above 90%. The 
imaging data saturate at bright magnitudes, so we restrict 
to sources with R > 20, leaving 161 photometric redshifts. 
The photometric redshifts {open triangles) for sources that 
do not also have spectroscopic redshifts are shown in Fig- 
ure 1, as is the completeness of our HDF-N combined 
photometric and spectroscopic redshift sample versus R 
magnitude {dashed histogram). Eighty-three sources with 
i? > 20 have both spectroscopic and photometric redshifts, 
and 92% of these agree to within a multiplicative factor of 
1.25. This agreement is independent of spectral type, ex- 
cept for the very small number of broad-line AGNs, where 
the failure rate is higher. The success rate of the photo- 
metric redshifts is consistent with that expected from the 
assigned odds. Since we do not have photometric redshifts 
for the SSA13 sources, and since photometric redshifts do 
not contain spectral type information, we only use the pho- 
tometric redshifts to test the effects of completeness on the 
total LF. 

For the spectroscopically identified sources, we classified 
the optical spectra into four spectral types, roughly fol- 
lowing the procedure used by Sadler et al. (2002, hereafter 
S02) to analyze low redshift 1.4 GHz samples. We clas- 
sified sources without strong emission lines (EW([OII])< 
3 A or EW(Ha-l-NII) < 10 A) as absorbers; sources with 
strong Balmer lines and no broad or high-ionization lines 
as star formers; sources with [NcV] or CIV lines or strong 
[OIII] (EW([OIII] 5007 A) > 3 EW(H/3)) as Seyferts; and, 
finally, sources with FWHM> 1000 km as broad-line 
AGNs. The breakdown by spectral type in each field, and 
also by radio fiux, is given in Table 1. S02 argued that the 
absorbers must be AGN dominated since star formation 
cannot produce the synchrotron emission in these galax- 
ies. They combined all but the star formers into a single 
AGN class. 

The deep CDF-N X-ray data provide a check and a re- 
finement of the classifications. We computed the 2 — 8 keV 
rest- frame luminosity, L^, for each source with a spectro- 
scopic redshift. In Figure 2 we plot versus radio power 
for the AGNs (classes 3 and 4; triangles), the star form- 
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ers {diamonds), and the absorbers (squares). The sources 
classed as AGNs arc generally X-ray luminous. Of the 40 
AGNs, 29 are detected in the 2 - 8 keV band, 26 of which 
have > 10^^ ergs s~^. Such luminosities would securely 
identify them as AGNs (Zezas, Georgantopoulos, & Ward 
1998; Moran, Lehnert, & Helfand 1999). The remaining 
upper bounds are consistent with all of the sources having 
> 10^^ ergs s"!. 




1.4 GHZ LUMINOSITY (ergs s"' Hz"') 



Fig. 2. — 2 — 8 kcV luminosity versus 1.4 GHz power for the 
HDF-N radio sample. Triangles show sources with AGN optical 
spectra {large symbols — broad-line AGNs; small symbols — Seyferts); 
diamonds show star formers; squares show absorbers. Arrows de- 
note upper limits for sources not detected in the 2 Ms CDF-N. Solid 
diagonal line is an approximate linear fit to the AGN population 
and is just shown for guidance. Dashed horizontal line shows the 
luminosity that would securely identify a source as an AGN. 

By contrast, nearly all of the star formers have low X-ray 
luminosities. Only six of the 55 sources are detected in the 
2 — 8 keV band, and only four have > 10*^ ergs s~^. 
The absorbers are the most unusual class. If they are 
AGN dominated, then they are very underluminous in 
their rest-frame 2 — 8 keV luminosities relative to their 
radio power. This suggests that they could be Compton- 
thick or near Compton-thick sources. It may therefore 
be useful to maintain the distinction between the AGNs 
with clear optical signatures and the absorbers. However, 
to be consistent with S02, we combine the two classes. 
The absorbers class has about half as many sources as the 
spectroscopic AGN class and a similar distribution of radio 
power. Thus, the AGN radio LF may be roughly corrected 
to a spectroscopic AGN LF by multiplying by 0.66. Some 
sources that have star- forming optical signatures may also 
correspond to this type of source: X-ray quiet, radio de- 
tected AGN with no spectroscopic AGN signatures. There 
may be ongoing star formation in the galaxy, but the AGN 
dominates the radio power. Thus, the optically classified 
star formers may only be an upper limit to the population 
where star formation dominates the radio power. 

For our final classification, we classed all absorbers, 
star formers known to have measured (not limits) > 
10^^ ergs s~^, and sources with AGN optical spectra as 
AGNs. We classed the remaining sources as star formers. 
This enables straightforward comparisons with the local 
results of Condon (1989), Condon, Cotton, & Broderick 
(2002), Machalski & Godlowski (2001, hereafter MGOl), 
and S02. 



3. DISCUSSION 

We constructed the LFs for the 1.4 GHz power over the 
range z = 0.5 — 1.5 using the 1/Va method (e.g.. Felten 
1977). This is shown for the total spectroscopic sample 
(AGNs and star formers; triangles) in Figure 3a, where it 
is compared with the local measurements of Condon (1989; 
solid line). To investigate the effects of incompleteness, 
we computed the maximum total LF (squares) for just the 
HDF-N using the much more complete combined photo- 
metric and spectroscopic samples and placing all sources 
without a photometric or a spectroscopic redshift into the 
z = 0.5 — 1.5 bin. Many of these probably lie at higher 
redshifts, so this is an upper bound. Finally, we assumed 
that the HDF-N sample was only 50% complete in the 
40 — 80 /LtJy range, which approximates the completeness 
corrections computed by Richards (2000). The maximum 
correction for incompleteness is a factor of 2.6 at the lower 
radio powers but is relatively small at the higher radio 
powers. 

In order to parameterize the evolution, we assumed that 
the local portion of the LF allocated to star formers in 
Condon (1989) undergoes pure luminosity evolution and 
that allocated to AGNs undergoes pure number density 
evolution. We used these assumptions to calculate the ex- 
pected LF in the 2 = 0.5 — 1.5 interval. (The use of the S02 
or the MGOl measurements, which were constructed us- 
ing more similar methodologies to ours for classifying the 
galaxies, does not change the diagrams significantly.) This 
simple model provides a reasonable fit to our measured LF, 
with the AGN number density evolving as (1-1-^)^, and the 
star-forming luminosities evolving as either (1 + z)^ (dotted 
line), if we match the observed LF, or as (1 + z)'^-^ (dashed 
line), if we match the maximum total LF. This systematic 
uncertainty is much larger than the statistical errors and 
primarily affects the evolution of the star-forming LF be- 
cause the uncertainties are larger at the lower luminosities 
where this term dominates. Our results for the evolution of 
the star-forming LF are consistent with those determined 
from local samples by MGOl and Condon et al. (2002), 
who found a similar evolution of (1 4- z)^^^ . Our results 
are also broadly consistent with those of HOO, who found 
a best pure luminosity evolution model with (1 -|- z)^'^. 

In Figures 3b and 3c we show the measured LFs in two 
narrower intervals (z = 0.6 — 0.9 and z ~ 0.9 — 1.4, re- 
spectively, chosen to match HOO) and by spectroscopic 
class (diamonds — star formers; squares — AGNs). For the 
star formers, we also show the results of HOO (triangles), 
which, within the uncertainties, agree with the present 
measurements. Our results arc again broadly consistent 
with the evolution described above (see Fig. 3b caption 
for description), although the AGN LF is slightly high in 
the z = 0.6 — 0.9 interval. This suggests that the power 
law evolution model may be too simple and that the evo- 
lution of the AGN LF is faster at low redshifts and slower 
at high redshifts. The star formers begin to show a more 
extended high luminosity tail in the z = 0.9 — 1.4 interval. 
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Fig. 3. — (a) Radio LF at z = 0.5 — 1.5 for all the spcctroscop- 
ically identified radio sources (triangles), compared with the local 
determination from Condon (1989; solid line). Log denotes the base 
10 logarithm. The Icr uncertainties are Poissonian. Squares denote 
the maximum total LF. Vertical dotted line shows the equivalent 
radio power of an Lpiji = lO^'^ Lq source. Dotted curve shows the 
expected LF if the local star-forming LF undergoes pure luminosity 
evolution ((l + z)'^) and the local AGN LF undergoes pure number 
density evolution ((l + z)'^). Uncertainties are primarily systematic, 
and luminosities evolving as (1 + z)'^'* provide a better fit {dashed 
line) to the maximum total LF. (b) Radio LF at 2 = 0.6 — 0.9 for the 
star formers {diamonds) and the AGNs {squares), to be compared 
with the local star former {steep solid curve) and AGN {shallow 
solid curve) determinations from Condon (1989). Triangles show 
the measurements of HOG (with a small correction to make the ge- 
ometries consistent) for the star formers. Dot-dashed curve shows 
the expected LF if the local star-forming LF undergoes pure lu- 
minosity evolution ((1 + z)^-*). Dashed curve shows the expected 
LF if the local AGN LF undergoes pure number density evolution 
((1 + z)3). (c) As for (b), but for the z = 0.9 - 1.4 redshift interval. 



The rapid evolution of the radio LF impHes that tlic 
number density of very luminous sonrc;es rises rapidly in 
the z = — 1 interval for both the star former and AGN 
classes though it may be begining to slow at redshifts 
z > 1. In Figure 4 we show the number density of star 
formers with ULIG radio power (solid squares) in three 
redshift intervals. The value from the local star-forming 
radio LF is given by the open square. The evolution is 
reasonably well described by a (1 + z)'' evolution {dashed 
line). This evolution is similar to that inferred by Barger 
et al. (2000) from a comparison of the local number density 
of ULIGs and near-ULIGs to the z = 1 — 3 number density 
of > 6 mJy radio/submillimeter sources (solid circle). If 
most of the submillimeter sources are dominated by star 



formation (Barger et al. 2000; Chapman et al. 2003), the 
match to the radio evolution is reasonable, if the radio 
galaxies classified as star formers are indeed dominated by 
star formation and the FIR-radio correlation holds to these 
redshifts as suggested by recent results (Garrett 2002). 
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Fig. 4. — Number density of star formers with ULIG radio power 

in the redshift intervals z = 0.1 — 0.6, z = 0.6 — 1.0, and z = 1.0 — 1.4 
{solid squares). Open square denotes the value from the local star- 
forming radio LF. Solid circles denote the local number density of 
ULIGs and near-ULIGs and the z = 1—3 number density of > 6 mJy 
radio/submillimeter sources (from Barger et al. 2000). Dashed line 
shows a (1 -|- z)'' evolution. 



4. SUMMARY 

We have given a precise determination of the high power 
end of the radio LF at high redshifts (z ^ 1). The LF was 
shown to evolve rapidly with redshift, both for galaxies 
with AGN spectra and for those with only star-forming 
signatures. This result is consistent with model expecta- 
tions based on the local LF and the radio number counts 
(e.g., Condon 1989). The number of sources with radio lu- 
minosities that would correspond to ULIGs (based on local 
normalizations) matches the observed evolution from the 
local ULIG population to the radio/submillimeter sources 
at z > 1. 
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from NASA grant HST-GO-09425.30-A, NSF grant AST- 
0084847, NSF CAREER award AST-0239425, the UW Re- 
search Committee, the Alfred P. Sloan Foundation, and 
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Table 1 

Breakdown of spectral type by field and radio flux ((uJy) 



CATEGORY 


HDF 


SSA13 


ALL 


25-50 


50-100 100-200 > 200 


total 


195 


151 


346 


37 


186 


77 


46 


identified 


115 


82 


197 


20 


98 


51 


28 


stars 





2 


2 





2 








absorbers 


15 


14 


29 


1 


16 


4 


8 


star formers 


57 


49 


106 


10 


51 


33 


12 


Seyferts 


36 


13 


49 


8 


23 


12 


6 


broad-line 


5 


3 


8 





5 


1 


2 



